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We present the results of a search for B + — > 7^ + ^, where £ = e, fi. We use a sample of 232 
million BB pairs recorded at the T(45) with the BABAR detector at the PEP-II B Factory. We 
measure a partial branching fraction AS in a restricted region of phase space that reduces the effect 
of theoretical uncertainties, requiring the lepton energy to be between 1.875 and 2.850 GeV, the 
photon energy to be between 0.45 and 2.35 GeV, and the cosine of the angle between the lepton and 
photon momenta to be less than —0.36, with all quantities computed in the T(4S) center-of-mass 
frame. We find AB(B+ -> -y£+v e ) = (-0.3t^(stat) ± 0.6(syst) ± O.l(th)) x 10~ 6 , assuming lepton 
universality. Interpreted as a 90% C.L. Bayesian upper limit, the result corresponds to 1.7 x 10~ 6 
for a prior flat in amplitude, and 2.3 x 10 -6 for a prior flat in branching fraction. 

PACS numbers: 13.20.He, 13.30.Ce, 12.38.Qk, 14.40.Nd 



At tree level, the branching fraction (BF) for radiative 
leptonic B decays is given by: 



B{B+ -► -ftvi) 



G F \Vub\ f 2 5 

288tt 2 JB B 



As 



nib 



(1) 

where tub is the B + meson mass, mi, is the MS b quark 
mass, tb is the B + meson lifetime, fs is the B meson 
decay constant, Qi is the charge of quark flavor i and Xb 
is the first inverse moment of the B light-cone distribu- 
tion amplitude [3, H, a quantity that enters into theo- 
retical calculations [3J] of the BF of hadronic B decays 
such as B — > 7T7T, and is typically taken to be of the order 
of Aqcd- Thus, a measurement of B(B + — ► ^l^vi) can 
provide a determination of As free of hadronic final-state 
uncertainties. The best current 90% C.L. upper limit on 
the full BF is 5.2 x lCT 5 [il, for B+ -> 7/x+^. 

However, Eq.(l) is based on the assumption that the 
factorization relation for the vector and axial- vector form 
factors is valid over the entire phase space. Instead, one 



samples, to discriminate signal from background. We 
then extract the number of signal events in data using a 
binned maximum-likelihood (ML) fit. 

The backgrounds are divided into three categories: 
continuum (non- £?£?), specific exclusive b — ► ulvg decays, 
and "generic B" decays, defined as a combination of all 
B hadronic decays, b — > c£vi decays, and the remain- 
ing inclusive b — > utvi decays. In particular, we study 
the seven exclusive b — > u£vg modes: B + — > h°£ + ve 
(h° = 7r°,p°,77,7/,w) and B° -> h~l + v t (h~ = n~,p-), 
referred to below, for each h, as the u h mode" . 

Our signal MC samples were generated using the tree- 
level model of Ref. [l|. The 7r° and 7r ± mode samples 
were generated using the form factor parameterization of 
Ref. [8j , with the value of the shape parameter based on 
lattice QCD results @. Light cone sum rule-based form 
factor models were used to generate samples for the p° 
P 



^ , and uj modes [lfj , and r\ and 77' modes 11 1 
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We find an excess of events in the off-peak data com- 

► qq,r + T~, and in the 



can relate, at tree-level, X B to a_partial BF, A/3, over a P ared to continuum MC (e+e 
restricted region of phase space 



AB 



G 2 F \V ub \ 

1 327T 4 



■f B T B m 3 B [a + bL + cL 2 ], (2) 



where L = (7Bb/3)(1/Ab + l/(2m(,)), the first term de- 
scribes the effects of photon radiation from the lepton, 
the third term the internal photon emission, and second 
their interference. The constants a, b, and c can be pre- 
dicted model-independently using factorization at large 
photon energy, the kinematic region for our analysis. 

We present herein the results of a search for charged 
B meson decays B + — ► j^vi, where I = e,/x ("elec- 
tron channel" , "muon channel" )Q. Our measurements 
are based on a sample of 232 million BB pairs recorded 
with the BABAR detector Q at the PEP-II asymmetric- 



muon channel, /z + /Lt~7), with the excess more pronounced 
in the electron channel. This is likely to result from 
unmodclcd higher-order QED and hadronic two-photon 
events. We thus use off-peak data instead of continuum 
MC to represent continuum background in our analysis. 

We take as the signal lepton and photon the high- 
est center-of-mass (CM) energy electron (muon) and the 
highest CM energy photon candidate in each event. The 
remaining charged tracks, each assigned a pion mass, and 
neutral clusters, treated as photons, are assigned to the 
"recoil B" candidate. We reconstruct the recoil B in two 
ways: we construct an "unsealed" recoil momentum as 
the sum of the CM 3-momenta of its constituents, and 
we define a "scaled" recoil momentum in the direction of 
the unsealed recoil, with its magnitude determined from 



energy e+e" storage rings, comprising an integrated lu- the CM energy of the T(AS) and the B ± mass. Using 
minosity of 210.5 fb _1 collected at the T(AS) resonance either the scaled or unsealed momentum, we reconstruct 
("on-peak"). We also use 21.6 fb" 1 recorded approxi- th c 3-momentum of a corresponding scaled or unsealed 

signal neutrino candidate. The reconstructed neutrino 
CM energy is calculated as the difference between the 
CM beam energy and the sum of the lepton and photon 
candidates' CM energies. 

We optimize a set of selection criteria for the best sig- 
nal sensitivity at a significance of 3cr using MC samples, 



mately 40 MeV below the T(45) ("off-peak"). 

The analysis procedure consists of selecting a lepton 
and photon recoiling against a reconstructed B, and iden- 
tifying signal candidates by reconstructing the neutrino 
using missing energy and momentum. We use a variety 
of selection criteria, optimized using Monte Carlo (MC) 
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splitting each in half, with one sample used for the opti- 
mization and the other used to evaluate its performance. 

On the signal side, we require that the electron (muon) 
have a CM energy between 2.00 and 2.85 (1.875 and 
2.775) GeV. We require that the photon have a CM en- 
ergy between 0.65 and 2.35 (0.45 and 2.35) GeV. We 
define cos6^ 7 to be the cosine of the angle between 
the lepton and photon in the CM frame, and require 
its value to be less than —0.42 (—0.36). We require 
-1.10(-1.05) < cos9by < 1.10(1.00), where cos6 B y 
is the cosine of the angle between the signal B and the 
lepton-photon combination Y in the CM frame [13] , com- 
puted from the known B mass, the beam energy, and the 
3-momenta of the signal lepton and photon. 

In order to reduce background from neutral hadrons, 
we require the lateral moment [l3T ] of the electromagnetic 
calorimeter energy distribution of the signal photon can- 
didate to be less than 0.55 for both channels. The polar 
angle of the photon candidate in the laboratory frame 
is required to be between 0.326 and 2.443 rad for both 
channels. We pair the candidate with every other neutral 
cluster in the event and reject events with a pair invariant 
mass in the tt° mass range 123-147 (116-148) MeV. 

We require the difference between the total CM energy 
of the recoil B constituents and the CM beam energy to 
be between -5.0 and 0.9 (-2.5 and 0.7) GeV. For the 
neutrino reconstruction, we require that both the scaled 
and unsealed neutrino polar angle in the laboratory frame 
be between 0.300 and 2.443 rad for both channels. 

To reduce continuum background, we require the ratio 
of the second to zcroth Fox- Wolfram moment 14| of all 
charged tracks and neutral clusters to be less than 0.5, 
and the absolute value of the cosine of the angle between 
the CM thrust axes of the recoil B and the lepton-photon 
system be less than 0.98 (0.86). We use a Fisher discrim- 
inant, T = 00^0 + 02^2, calculated from the momentum- 
weighted zeroth and second Legendre moments, Lq and 
L-2 , of the recoil B about the lepton-photon CM thrust 
axis, with coefficients ao and ai equal to 0.43 and —1.86 
(0.008 and —1.590), respectively T is required to be 
greater than 1.50 (0.310). 

In the electron channel, we veto two-photon events 
via the charge-angle correlation of the signal lepton aris- 
ing from the initial state. For a positively (negatively)- 
charged signal electron, we require the cosine of its CM 
polar angle to be between —0.74 and 0.78 (—0.94 and 
0.70). In the muon channel, we require this variable to 
be between —1.00 and 0.78 for both charges. These crite- 
ria were optimized on a loosely-selected sample of events, 
where the off-peak data are used for the continuum, and 
the MC for the signal and other backgrounds. 

We also reject two-photon events using a parameter- 
ized combination of the missing CM momentum in the 
beam direction and the invariant mass of the hypotheti- 
cal two-photon system. For the muon channel, the entire 
observed event is taken as the two-photon system, while 
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FIG. 1: Electron-channel Aep vs. mES signal MC, using a 
color scale to represent relative contents of each bin. 

for the electron channel, the signal electron is assumed 
to be from the initial state, and so is excluded from the 
two-photon system. The selection criterion was adjusted 
to preserve a 94% efficiency for signal for both channels. 

After applying our selection criteria, we use the two- 
dimensional distribution of Aep, the difference between 
the scaled neutrino candidate's CM energy and the mag- 
nitude of its 3- momentum, and toes, the invariant mass 
of the recoil B, calculated from its unsealed CM 3- 
momentum and the CM beam energy, as inputs to the 
ML fit. These distributions provide distinct signatures 
for signal, B background, and continuum, with the sig- 
nal distribution shown in Fig.[TJ The signal (S) and three 
sideband (Bl, B2, B3) regions were selected to maximize 
separation of signal from BB and continuum background. 

We extract signal events by fitting on-peak data for the 
contributions of signal and background, while allowing 
the predicted shapes of signal and background to vary 
within statistical uncertainties. The scale of signal and 
generic B contributions are allowed to vary, while the 
scale of off-peak data is fixed using the on-peak/off-peak 
luminosity ratio. For the seven scmilcptonic (SL) modes, 
we fit for three of the BFs and relate the other four to 
them as follows: The 7r ± and p^ mode BFs are obtained 
from BaBAR measurements 12|, and the r\ mode BF is 
obtained from CLEO [l5|. The charged and neutral 7r 
and p modes are related by the lifetime ratio, t b ± /t b o = 
1.071 ± 0.009 [3, and an isospin factor of 2. The uj 
mode BF is taken as equal to the p° mode BF. We take 
the ratio of the 77 to rf mode BFs to be 2.057±0.020[l7|. 

We maximize a likelihood function consisting of the 
product of four Poisson probability distribution functions 
(PDFs), modeling the total counts in each of the four 
regions, three Gaussian PDFs for the BFs of the three 
SL modes, and 40 Poisson PDFs for the 4-rcgion shapes 
of the various samples. All of the shapes are obtained 
from MC, except for continuum, where off-peak data are 
used, introducing a larger statistical uncertainty. 

Each Poisson PDF that models the total count in one 
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of the four fit regions has a measured value obtained from 
the on-peak data count, and an expected value based 
on the fitted contributions of signal and background, in- 
cluding fitted variations of the shapes. For the seven 
SL modes (where three of the fitted BFs are indepen- 
dent), the variances of the three Gaussian likelihoods are 
obtained from the published statistical and experimental 
systematic uncertainties, combined in quadrature. In all, 
there are 47 PDFs, and 45 free parameters. 

We fit for the partial BF AB for the kinematic region 
with lepton CM energy between 1.875 and 2.850 GeV, 
photon CM energy between 0.45 and 2.35 GeV, and 
cos#£ 7 less than —0.36 — the union of the electron and 
muon channel regions. We perform three fits: separate 
electron and muon channel fits, and a joint fit in which 
the signal and three SL BFs are constrained to be equal 
for the two channels. For each fit, errors on the fitted sig- 
nal BF are obtained by finding the two values at which 
the signal BF likelihood decreased by a factor of e -1 / 2 . 
Table U shows the results from the joint fit. 

TABLE I: Comparison of fit results and experimental obser- 
vations for the joint fit to the muon and electron channels. 
For each of the four fit regions, the individual fitted contri- 
butions from continuum (cont.), BB background, and signal 
are shown, along with their total. The on-peak and off-peak 
(scaled to the integrated on-peak luminosity) observations are 
shown for comparison, in indented rows, and are not included 
in the "Total fit" value shown. 

Muon channel 

S Bl B2 B3~ 

Fit cont. 20.0T11.8 116.3T14.7 42.6T12.8 213.2T42.1 

Off-peak 23.0T16.2 158.1T40.8 17.4T12.3 219.7T45.8 

Fit BB 59. 1± 8.5 61. 0± 9.9 61. 7± 9.8 286.6T46.6 

Fit signal -5.2T13.8 -1.3± 3.4 -0.4T 1.0 -0.2± 0.5 

Total fit 74.0± 8.1 176.0T12.4 103.9T 9.8 500.0T22.1 

On-peak 73.0T 8.5 170.0±13.0 111.0±10.5 498.0T22.3 

Electron channel 

S Bl B2 B3~ 

Fit cont. 55.4T20.5 181.1T16.2 48.9T14.1 356.7T54.4 

Off-peak 41.4T20.7 239.7T48.9 79.0T27.9 294.5T52.9 

Fit BB 69. 2± 8.5 59.2T 8.5 140.1T15.5 393.8T57.2 

Fit signal -8.4T22.3 -1.5± 3.9 -1.2± 3.3 -0.4± 1.0 

Total fit 116.2T10.3 238.7T14.5 187.7±12.5 750.2T26.5 

On-peak 119.0T10.9 231.0T15.2 176.0T13.3 764.0T27.6 



Table [XT] shows all systematic uncertainties on AB ex- 
cept for theoretical uncertainties on the signal model, 
which are shown in Table Hill 

The experimental systematic errors result from un- 
certainties on the data/MC consistency with respect to 
tracking efficiency of the signal lepton, particle identifica- 
tion efficiency of the signal lepton, reconstruction of the 
signal photon energy, selection criteria efficiency uncer- 
tainties, and uncertainties on the data/MC consistency 



in the shape of Aep and toes ■ All of these were evaluated 
using a number of control samples, including /i + /i~7, 
e+e-7, and B+ -> 7r + U°(-> K+n-). 

The uncertainty on the number of produced BB pairs 
is 1.1%. In our fits, we further assumed a charged-to- 
neutral B production ratio of 1.0, and determine the 
systematic uncertainty by varying within the measured 
interval 1.020 ± 0.034 [l|. 

For systematic errors due to the theoretical uncertain- 
ties on the 71" and p mode BFs and form factor models, 
we refit, applying correlated variations in the BFs and 
form factor models, and take the magnitude of the largest 
change in signal BF as the associated systematic. 

For systematic errors arising from the theoretical un- 
certainties on the r\ and rj' mode BFs, we vary the as- 
sumed BFs by ±10%[l5j and refit to obtain a systematic 
error. We find a negligible systematic from the uncer- 
tainty on the ratio of the r\ to rj' mode BFs. 

In the generic B sample, a significant fraction of events 
contributing to our fit are non-resonant B — > X u £ + i/£ 
events. We obtain the systematic error due to uncer- 
tainty on the B — » X u t + vi BF by fixing the total contri- 
bution of generic B decays, as predicted by MC. 

TABLE II: Systematic uncertainties on AS. All additive sys- 
tematic values have been multiplied by 10 6 . 



Multiplicative 


Muon Electron Joint 


Tracking efficiency 


1.3% 


1.3% 


1.3% 


Particle ID 


3.5% 


2.2% 


2.1% 


Neutral reconstruction 


1.6% 


1.6% 


1.6% 


Selection efficiency 


6.0% 


5.0% 


6.0% 


B counting 


1.1% 


1.1% 


1.1% 


Charged to neutral B ratio 


9.4% 


9.4% 


9.4% 


Additive 


Shape of Aep vs. tties 


0.3 


0.2 


0.3 


77 mode BF 


0.3 


0.1 


0.2 


7r,p mode BF, ff 


0.3 


0.4 


0.1 


B X u l+v t BF 


0.4 


0.2 


0.3 



The theoretical uncertainty within the kinematic re- 
gion of AB is conservatively estimated by evaluating the 
change in efficiency when the model of Ref. [l[ is modified 
by setting the axial vector form factor equal to zero. 

The results for AB are given in Tables UTTl and HVl We 
determine 90% C.L. Bayesian upper limits by integrating 
the signal BF likelihood with two different priors, both 
of which take values of for negative values of the signal 
BF: a prior flat in the BF ("flat BF prior"), and a prior 
fiat in the square root of the BF ("flat amplitude prior"), 
equivalent to assuming a flat prior for \V u b\ or fp- 

For our kinematic region, the constants a, b, and c 
of Eq.(2) are 0.88, -3.24, and 3.25, respectively @. Us- 
ing input values of f B = 216MeV [lj], \V ub \ = 4.31 x 
10~ 3 [16], t b = 1.638 ps [lj], and m b = 4.20GeV 0, 
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TABLE III: Comparison of AB two-sided results for all three 
fits. All values have been multiplied by 10 6 . 

All Central Statistical Systematic Theoretical 

values xlO 6 value uncertainty uncertainty uncertainty 



EIF (European Union) and the A. P. Sloan Foundation. 



Muon 

Electron 

Joint 



-1.33 
0.11 
-0.25 



+1.74 
-2.20 
+1.73 
-2.13 
+1.33 
-1.53 



+0.80 
-0.87 
+0.61 
-0.59 
+0.60 
-0.64 



0.03 
0.08 
0.07 



TABLE IV: The 90% Bayesian upper-limits for all three fits, 
for the two different choices of prior, in terms of AB. 

Prior flat in amplitude Prior flat in BF 



Muon 

Electron 

Joint 



< 1.5 x 10 

< 2.2 x 10 

< 1.7 x 10 



< 2.1 x 10" 

< 2.8 x 10" 

< 2.3 x 10" 



our 90% C.L. Bayesian limits on AB correspond to val- 
ues of A B of > 669 MeV and > 591 MeV, for the choice 
of the flat amplitude and flat BF priors, respectively. 

Given a theoretical model, a measurement of AB may 
be converted into an estimate of the total BF. In the 
model of Ref. [![, the result of the joint fit corresponds 
to a BF of (-0.6t^(stat)±^(syst)) x 1CT 6 , and 90% 
C.L. Bayesian upper limits of 3.8 x 10~ 6 and 5.0 x 10~ 6 
for the flat amplitude and flat BF priors, respectively. 
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